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Abstract 



The stability of equilibrium configurations of galaxies or stars are time 
honored problems in astrophysics. We present mathematical results on 
these problems which have in recent years been obtained by Yan Guo and 
the author in the context of the Vlasov-Poisson and the Euler-Poisson 
in - model. Based on a careful analysis of the minimization properties of 

' _ ^ , conserved quantities — the total energy and so-called Casimir functionals — 

1 non-linear stability results are obtained for a wide class of equilibria. 

o _ 

^ ■ 1 Introduction 

q . Under suitable idealizing assumptions the time evolution of a galaxy can be 

[ modeled by the Vlasov-Poisson system 

CC 

d t f + vV x f-VU-V v .f = 0, 
AU = 4-irp, lim U(t,x)=0, 

. W ; P(t,x)= l f{t,x,v)dv. 



Here f = f(t,x,v)>0 denotes the density of the stars in phase space, t&WL, x,v S 
1R 3 stand for time, position, and velocity, p = p(t,x) is the induced spatial mass 
density, and U — U(t,x) is the gravitational potential of the galaxy. 

The problem we wish to address is the non-linear stability of stationary 
solutions of this system. Our approach will automatically address this question 
for the Euler-Poisson system as well, which describes a self-gravitating fluid ball, 
i.e., a barotropic star. The latter model is presented in Section 4. 

By definition, a given steady state /o is stable if for any neighborhood TV of fo 
there exists another neighborhood M of fo such that any solution of the system 
starting in M will remain in N for all times. This is the usual mathematical 
definition of Lyapunov stability. In the case of an infinite dimensional dynamical 
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system such as the Vlasov-Poisson system the choice of the proper concept of 
"neighborhood" is a non-trivial part of the stability problem. We emphasize 
that no approach of the solution to the steady state is asserted — that would 
be the concept of asymptotic stability. Since the system as stated does not 
include dissipative effects an approach to a particular steady state is not to be 
expected in a strict sense; we do in the present notes not enter into the highly 
interesting questions of course graining, phase mixing, etc. The existence of 
global-in-time solutions of the system under consideration at least for initial 
data close to the steady state is an integral part of the stability assertion. For 
the Vlasov-Poisson system it was shown in pQ that every smooth, compactly 
supported initial datum for / launches a unique global-in-time smooth solution; 
a fairly short proof due to J. Schaeffer is given in [2]. 

There is a vast astrophysics literature on the stability question. However, 
essentially all investigations that we are aware of proceed via linearization. This 
approach suffers from at least two major difficulties: Firstly, there is no gen- 
eral theory for infinite dimensional dynamical systems as to how to pass from 
linearized to non-linear stability. Secondly, it is well known that if A is an eigen- 
value, i.e., the linearized system has a solution of the form e xt g(x,v), then the 
same is true for —A. Hence the optimal case regarding stability occurs if all such 
eigenvalues are purely imaginary, which is precisely the case where stability for 
the non-linear system does not follow, not even in finitely many dimensions. 

We will prove non- linear stability of certain steady states by identifying them 
as minimizers of a conserved quantity in terms of which the above neighborhoods 
N and M are then defined. More precisely, for a state / = f(x,v) > we denote 
the induced spatial mass density and potential by 

Pf {x):= J f(x,v)dv, U f (x):=-J 



#kin(/) := I I I \v\ 2 f{x,v)dvdx, (1.1) 



and we introduce the functionals 
1 
2 

E P at(f) :=\j U f (x)pf(x)dx = ~ J \VU f (x)\ 2 dx, (1.2) 

the kinetic and the potential energy of the state /. The total energy 

H:= Eki n + i?p t 

is conserved along solutions of the Vlasov-Poisson system, but it is indefinite, 
and it has no critical points, i.e., the linear part in an expansion about any state 
/o with potential Uq does not vanish: 

H(f) = H(f ) + J I Qm 2 + CA)) {f-fo)dvdx--^ J \VU f -VU \ 2 dx. 

However, according to Liouville's theorem the characteristic flow correspond- 
ing to the Vlasov equation preserves phase space volume, and hence for any 
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reasonable function 4> the so-called Casimir functional 

C(f):= JJ<S>(f(x,v))dvdx 

is conserved as well. If we expand the energy-Casimir functional 

H C :=H + C 

about an isotropic steady state 

fo(x,v) = 4>(E), E = E(x,v):=^\v\ 2 + U (x), 

we find that 

H c (f)=n c (fo)+ ff(E+*Va))tf-fo)dvdx 

/ \VUf-VU \ 3 dx+^ J I &'(f )(f-fo) 2 dvdx+.... 

At least formally, we can choose $ such that /o is a critical point of He, namely 
$' = — provided <j> is invertible. The essential problem now is the following: 
In order for the steady state to have finite total mass the function (j> must vanish 
above a certain cut-off energy. For <\r Y to exist cf> should thus be decreasing, 
at least on its support. But then $" is positive and the quadratic part in the 
expansion indefinite. Since one would like to use this quadratic part for defining 
the concept of distance or neighborhood, the method seems to fail. This state 
of affairs had been observed by various authors, with the proposed conclusion 
that the energy-Casimir method does not work for the stellar dynamics case of 
the Vlasov-Poisson system, cf. for example [3]. If the issue is the stability of a 
plasma, the sign of the source term in the Poisson equation and hence also the 
one in front of the potential energy difference in the expansion above is reversed, 
and up to some technicalities stability follows, cf. 0]. 

The approach developed by Yan Guo and the author to overcome this diffi- 
culty is as follows. Starting with a given function $ which defines the Casimir 
functional we try to minimize the energy-Casimir functional Tic under the con- 
straint that only states with a prescribed total mass M > are considered. Un- 
der suitable assumptions on $ a minimizer does exist in spite of the fact that 
the quadratic term in the expansion above is indefinite. One can then show that 
such a minimizer is a non-linearly stable steady state. The exact statements of 
these results will be given in the next section — the main assumption is that $ is 
strictly convex which is equivalent to fo(x,v) = <fi(E) being strictly decreasing. 

The crucial step is to prove the existence of a minimizer. Here we first 
construct out of the energy-Casimir functional Tic a reduced functional Ti r 
which is defined on the space of spatial mass densities p in such a way that there 
is a one-to-one correspondence between the minimizers of the two functionals. 
This reduced functional is analyzed in Section 3. The original motivation for 
introducing it was purely mathematical: It is defined on a simpler space, and 
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the troublesome part of the original functional is the quadratic and negative 
definite potential energy, but the latter depends on the spatial mass density pf 
and not directly on /. The detour via the reduced functional has a beautiful 
pay-off: The minimizers of the reduced functional are stable steady states of 
the Euler-Poisson system with a macroscopic equation of state corresponding to 
the microscopic equation of state f Q = <j>(E) induced by the Casimir functional. 
Hence via this reduction procedure we obtain a non-linear proof of what is 
often referred to as Antonov's First Law 5, p. 305]: A spherical stellar system 
with fo = fo(E) and dfo/dE<0 is stable if the barotropic star with the same 
equilibrium density distribution is stable. This relation to the Euler-Poisson 
system, i.e., to the stability of gaseous stars, is investigated in Section 4. 

In Section 5 we give the main arguments leading to the existence of a min- 
imizer for the reduced functional. Mathematically, this is the essential and 
non-trivial part; it can be skipped without compromising the understanding of 
the rest. 

To keep the presentation reasonably simple we restrict ourselves mostly to 
spherically symmetric, isotropic steady states. However, the method has also 
been applied to non-isotropic steady states, to axially symmetric ones, and to 
disk-like steady states. Some comments on these extensions together with other 
related results as well as open questions are collected in the last section. 

To conclude this introduction we should mention that none of the results 
presented here are new, although the way they are presented is new. The 
motivation for these notes is to collect in one place the main features of our 
method, and to present them in such a way that the readers can hopefully 
appreciate the ideas involved. For some details which are not so relevant for 
the main argument we refer to existing papers, but mainly our presentation is 
aimed to be self-contained. We include almost no references to the astrophysics 
literature. This is really not done out of disrespect but due to the belief that our 
method is essentially the first to address the full non-linear stability problem 
for the systems under consideration. Should these notes inspire some comments 
or criticism from the astrophysics community we would truly appreciate this. 
Acknowledgments: These notes are an expanded version of my presentation 
at the workshop "Nonlinear Dynamics in Astronomy and Physics" in November 
2004 at the University of Florida. I truly appreciated the kind invitation to this 
inspiring event, as well as the feedback I received there. The results reported 
here originate from my collaboration with Y. Guo, Brown University, whom I 
would like to thank for many stimulating discussions. 

2 Nonlinear stability for the Vlasov-Poisson 
system — statement of the results 

We fix a Casimir functional C, i.e., a function $ such as 

$(/) = / 1+1/fe J>0, with 0</c< 3/2, (2.1) 
more generally: <& G C 1 ([0,oo[) with $(0) = = $'(0), and 
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(<I>1) $ is strictly convex, 

($2) $(/) > Cf 1+1 / k for / > large, with < k < 3/2, 

($3) $(/) < C*/ 1+1 / fc ' for / > small, with < fc' < 3/2. 

For a given constant M > we want to minimize the energy- Casimir functional 
Tic over the constraint set 

T M := {/ 6 Ll(R e ) | JJfdvdx = M, £y n (/) + C(/)<oo|. 

Here Li (M 6 ) denotes the set of non-negative, integrable functions on R 6 . One 
can show that the potential energy is defined on this set and the two forms of 
Spot given in l|1.2|) are equal, cf. 6 , Lemma 1] . Due to conservation of mass the 
constraint set Tm is invariant under solutions of the Vlasov-Poisson system. 
The main step in the stability analysis is to establish the following theorem: 



Theorem 1 The energy- Casimir functional Tic is bounded from below on Tm 
with Km '■= inf^M He < 0- Let (fj)cJ-M be a minimizing sequence of He, i-e-, 
7~Lc(fj) ~~ > h-M ■ Then there exists a function /oG-^Mj a subsequence, again de- 
noted by (fj) and a sequence (a,j) CK 3 of shift vectors such that for the induced 
gravitational fields 

T a >VU f] = Vt//. (- + aj ) - VU fo m i 2 (M 3 ), j - oo. 

The state /o minimizes the energy- Casimir functional: 1~Lc{f<y)—h'M- 

We shall see shortly that to conclude stability of the state /o the theorem 
is needed in the above form; mere existence of a minimizer is not sufficient. A 
proof via a reduced functional is given below. The main difficulty is seen from 
the following sketch of the argument: To obtain a lower bound for the functional 
on the constraint set is easy, and by Assumption ($2) minimizing sequences can 
be seen to be bounded in L 1+1/>k . A standard analysis result then implies that 
such a sequence has a weakly convergent subsequence, which means that for any 
test function g from the dual space L 1+k the convergence J fjg^> J fog holds, 
cf. Sec. 2.18]. The weak limit fo is the candidate for the minimizer, and one 
has to pass the limit into the energy-Casimir functional. This is easy for the 
kinetic energy, the latter being linear, and for the Casimir functional which is 
convex due to Assumption ($1) it relies on Mazur's lemma, cf. (7J Sec. 2.13]. 
The difficult part is the potential energy, for which one has to prove that the 
induced gravitational fields converge strongly in I? . This problem will we dealt 
with on the level of the reduced functional in Section 5. 

Since our minimization problem is invariant under spatial translations we 
obtain a trivial minimizing sequence by shifting a given minimizer in space. If 
for example we shift it off to spatial infinity we cannot obtain a subsequence 
which tends weakly to a minimizer, unless we move with the sequence. Hence the 
spatial shifts in the theorem arise from the physical properties of our problem. 
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The Euler-Lagrange identity corresponding to our constrained variational 
problem implies that any minimizer is a steady state of the Vlasov-Poisson 
system. For the proof we refer to |S] or [5]: 



Theorem 2 Let Jq^Tm be a minimizer with potential Uq. Then 




with Lagrange multiplier E$. In particular, fa is a steady state of the Vlasov- 
Poisson system. 

For example, the choice (|2.1|) yields the polytropic steady state 



It should be noted that the assumptions on <I> easily translate into assumptions 
on the steady state /o as a function of the particle energy, the main one being 
that this function is strictly decreasing on its support. Various additional prop- 
erties can be derived for these minimizers/steady states, in particular they are 
necessarily spherically symmetric, cf. [HI Thm. 3] or Thm. 2]. Non-symmetric 
steady states will be briefly considered in the last section. 

To deduce our stability result we expand Tic about the minimizer /q: 



with d(f,f ) = iff / = fo- This is due to the strict convexity of and the fact 
that on the support of f the bracket vanishes by Theorem |2 note also that 
IS if - fo) = for / e T M - For suitable $'s, we even have d(f, fo)>cff\f-fo\ 2 . 
The point now is that according to Theorem^the term with the negative sign in 
(|2.2|l tends to zero along any minimizing sequence. This allows us to use the sum 
of the two positive definite terms in the expansion as our measure of distance in 
the stability result. As mentioned in the introduction, initial data from the space 
Cc(M. 6 ) of continuously diffcrcntiable, compactly supported functions launch 
smooth global-in-time solutions of the Vlasov-Poisson system which preserve all 
the physically conserved quantities. As above, T a f(x,v) :— f(x + a,v). 

Theorem 3 For any e > there exists a 5 > such that for any solution ti—tf(t) 
of the Vlasov-Poisson system with f(0)GC^(M 6 )nJ r M the initial estimate 




Hc(f)-Hc(fo) = d(f,f )-— / \VU f -VU \ 2 dx (2.2) 




where for / £ Tm , 
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implies that for any t>0 there is a shift vector o£R 3 such that 

d(T a f(t)J ) + ±- J \T a WU m -VU \ 2 dx<e, t>0, 
(provided the minimizer /o is unique up to shifts). 

We will comment on the uniqueness assumption (and why it comes in paren- 
thesis) shortly, but first we give the proof of this result, which is surprisingly 
simple — the difficulty resides in the proof of Theorem 

Proof: Assume the assertion is false. Then there exist e > 0, tj > 0, fj(0) € 
C\ (R 6 ) n T M such that for j G N, 

d(M0)Jo) + ^ J \VU fj(0) -VU a \ 2 dx<± 

but for any shift vector a£l 3 , 

d(T a f 3 (t J )J ) + ±- J \T a VU fj{t]) -VU \ 2 dx>e. 

Since Tic is conserved, (|2.2[l and the assumption on the initial data implies that 
Hc(fj(tj)) =ftc{fj(0)) -»Wc(/o), i- c -> ifj( t j)) c ^ r M is a minimizing sequence. 
Hence by Theorem^ / |VJ7/.({.) — VUo\ 2 — >0 up to subsequences and shifts in x, 
provided that there is no other minimizer to which this sequence can converge. 
By Q2.2p . d(fj(tj),fo) — >0 as well, which is the desired contradiction. □ 
Some comments are in order: For the polytropic steady states one can show 
that for a given mass M the minimizer is indeed unique up to shifts, as assumed 
above, cf. Thm. 3]. In general, minimizers do not seem to be unique; for 
a numerically verified example of non-uniqueness we refer to Rem. 3 (b)]. 
However, minimizers always seem to be isolated up to shifts which is sufficient 
for the above statement to still hold true, cf. |UJ Thm. 3]. If there were a 
continuum of minimizers then this set of minimizers as a whole would be stable, 
cf. ^3 p. 242]. Finally, for a closely related approach to which we will come 
back in the last section it is shown in JT] that the assertion of Theorem 01 holds 
without /o being unique or isolated. We kept the former assumption to make 
the proof simple. 

The spatial shifts appearing in the stability statement are again due to the 
spatial invariance of the system: If we perturb /q by giving all the particles an 
additional, fixed velocity, then in space the corresponding solution travels off 
from /o at a linear rate in t, no matter how small the perturbation was. 

A nice feature of the result is that the same quantity is used to measure 
the deviation initially and at later times t. In infinite dimensional dynamical 
systems control in a strong norm initially can be necessary to gain control in 
a weaker norm at later times. Whether our concept of distance is appropriate 
from a physics point of view is open to debate — it is simply what comes out 
of the energy-Casimir method. For the polytropic steady states our approach 
has been extended to yield stability with respect to the L 1 -norm of /, cf. |12|. 
Definitely a weak point is the fact that the proof is not constructive: Given an 
e we do not know how small the corresponding S must be. 
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(v)dv = s> (3.1) 



3 The reduced variational problem 

We wish to factor out the dependence on the velocity variable in our minimiza- 
tion problem. Starting from a given function f = f{x,v) with induced spatial 
density pf = pf(x) we clearly decrease TLc{f) by minimizing for each point x 
over all functions g — g(v) which have as integral the value Pf{x). This procedure 
does not affect the potential energy and reduces the sum of the kinetic energy 
and the Casimir functional into a new functional which no longer depends on / 
directly but only on pf. More precisely, with 

*(*):= infjy Q| w | 2 ff + $( 5 )) dv \geL\iR 3 ), J g{ 

we have the estimate Hc(f) >W r (p/) where 

H r (p):= J V(p(x))dx + E pot (p). 

We now wish to minimize Ti. r over the constraint set 

ft M := jpe£+(M 3 ) | J*(p)<oa, Jp = My 

These constructions owe much to ^3]- Before we analyze the reduced problem 
we make sure that we can lift any information gained for the latter back to the 
level of the original problem: 

Theorem 4 For all JgFm the estimate Hc(f) >T~t r (pf) holds, with equality 
if f = fo is a minimizer. Let po £ TZm minimize Ti T and Uq = U Po . Then 

( (VTHEo-Uo) ,U <E , 
P °-{ ,U >E (3 - 2) 

with Lagrange multiplier Eq, and 

( (&)-HEo-e) ,e<e 0i i ||2xf , n 
fo-=y q \e>e . e=e ( x ' v > : = 2^ 

lies in Tm an d minimizes TLc ■ If on the other hand fo £ Tm minimizes Tic 
then pf £ TZm minimizes Tl r . 

Equation (|3.2|) is nothing but the Euler-Lagrange identity for the reduced prob- 
lem, cf. H|; the theorem is proven in detail in |14l Sec. 2]. 

Let us now consider the reduced variational problem in its own right. The 
function , J defining the reduced functional is taken from the following class: 
*£C 1 ([0,oo[), #(0) = = *'(0), and 
("fl) $ is strictly convex. 

(*2) #(p) > Cp 1+1 / n for p > large, with < n < 3, 
(*3) *(p) < Cp 1+1/n ' for p > small, with < n' < 3. 
In Section 5 we shall prove the following central result: 
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Theorem 5 The functional TL r is bounded from below on TZm. Let (pj) GTZm 
be a minimizing sequence ofTL r . Then there exists a sequence of shift vectors 
(ctj)clR 3 and a subsequence, again denoted by (pj), such that 



and po € TZm *s o, minimizer of TL r . 

We need to translate the conditions on \& back into conditions on <f>. To do so 
we denote the Legendre transform of a function h:M.— >] — 00,00] by 



This more explicit relation between $ and \& is established in [141 Sec. 2], and 
it allows us to show that the assumptions on $ imply the ones on Sfr if the 
parameters k and n are related by n = fc + 3/2, with the same relation holding 
for the primed parameters. 

Theorem 0] connects our two variational problems in the appropriate way 
to derive Theorem ^ from Theorem |SJ Firstly, He is bounded from below on 
Tm since this is true for TL r on TZm- Let (fj) cTm be a minimizing sequence 
for Tic- By Theorem H {p f] )c_TZ M is a minimizing sequence for TL r . Again 
by Theorem 0] we can lift the minimizer po of Ti r obtained in Theorem to a 
minimizer fa of Tic, and the proof of Theorem^is complete. 

Before we consider some of the ideas involved in the proof of Theorem we 
reinterpret it in terms of the Euler-Poisson system. 

4 Pay-off of reduction — The Euler-Poisson sys- 
tem 

If pqGTZm minimizes the reduced functional TL r , then po supplemented with 
the velocity field uq = is a steady state of the Euler-Poisson system 



T a3 P 3 ^Po weakly in L 1+1/ "(R 3 ), j^oo, 
T^V[/p 3 ^ V[/q strongly in L 2 (R 3 ), j^oo, 



h(X) :=sup refi (Ar-/i(r)). 



If \I/ arises from $ by reduction, i.e., by formula lj3.1Jl then 




d t p+V-(pu) = Q, 



pdtu + p(u-\7)u = —Vp~ pVU, 
AU = 4:iTp, lim U(t,x) = 0, 



with equation of state 



p = P(p):=p*'(p)-V(p). 
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This follows from the Euler-Lagrange identity 13.2fl . Here u and p denote the 
velocity field and the pressure of an ideal, compressible fluid with mass density 
p, and the fluid self- interacts via its induced gravitational potential U. This 
system is sometimes used as a simple model for a gaseous, barotropic star. The 
beautiful thing now is that obviously (po,uo = 0) minimizes the energy 

H(p,u):=- J \u\ 2 pdx + J $>(p)dx+E pot (p) 

of the system, which is a conserved quantity. Expanding as before we find that 

H(p,u)-H(po,0) = ^ J \u\ 2 pdx + d(p,p )--^ J \VU p ~VU a \ 2 dx 

where for p G TZm , 

d(p,p ):= ( [^(p)-^(po) + (Uo-E o )(p-po)}dx>0, 



with equality iff p = pq. The same proof as for the Vlasov-Poisson system implies 
a stability result for the Euler-Poisson system — the term with the unfavorable 
sign in the expansion again tends to zero along minimizing sequences, cf. The- 
orem 03 However, there is an important caveat: While for the Vlasov-Poisson 
system we have global-in-time solutions for sufficiently nice data, and these so- 
lutions really preserve all the conserved quantities, no such result is available 
for the Euler-Poisson system, and we only obtain a 

Conditional stability result: For every e > there exists a 6 > such that for 
every solution t<-^(p(t),u(t)) with p(0) £1Zm which preserves energy and mass 
the initial estimate 

l - J \u(0)\ 2 p(0)dx + d(p(0),p ) + ^ J \VU m -VU \ 2 dx<8 
implies that as long as the solution exists, 

i J \u(t)\ 2 p(t)dx + d{p(t),p ) + -±- J \VU p{t) -VU \ 2 dx<e 

up to shifts in x (provided the minimizer is unique up to such shifts). 

The same comments as on Theorem [21 apply also in this case. Because of 
the above caveat we prefer not to call this a theorem, although as far as the 
stability analysis itself is concerned it is perfectly rigorous. The open problem 
is whether a suitable concept of solution to the initial value problem exists. 

Now that minimizers of the reduced functional are identified as steady states 
of the Euler-Poisson system it is instructive to reconsider the reduction proce- 
dure leading from the kinetic to the fluid dynamics picture. First we recall that 
for the Legendre transform 

h'(0 = V KO + h( V ) =£r/^ {h)'(v) = t 
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If fo is a minimizer of Tic, 

Po = Jfadv = J (^'(Eo-Uo-^vAdv, 
P° = IJ \ v \ 2 fodv = J$(e q -U -^\v\ 2 ) dv. 
On the other hand, if po is a minimizer of 7i r , 

Pa = {V')- 1 (E -U ) = (W)'(E -Uo), 

Po = P(p )=p ^'(po) - *Oo) = *(* / (po)) =*(Eo - U ). 

In order for these relations on the kinetic and on the fluid level to fit we have 
to require that 



which is the relation between <I> and \& obtained by our reduction mechanism. 

5 Proof of the existence of minimizers for the 
reduced problem 

In this section we present the main arguments leading to the proof of Thcorcm|21 
Constants denoted by C may only depend on M and \& and may change their 
value from line to line. For a set Set 3 the indicator function 1<j equals 1 on S 
and vanishes outside. By Br we denote the ball of radius R about the origin. 
Step 1: Lower bound for H. r and weak convergence of minimizing sequences. We 
need to estimate the negative potential energy against the positive part of H r . 
The Hardy-Littlewood-Sobolev inequality 7, Thm. 4.3] tells us that 

-E pot (p)<C\\p\\l /5 . 

By the restriction on n, l<6/5<l + l/n, and by interpolation and (^2), 

/ f \ «/3 



Hence on TZm 

n/3 



n r (p)> J y{p)dx-C-c(^J q{p)dx^J 



(5.1) 



Since n<3 this implies that TC r is bounded from below on TZm' 

Km '■— inf Ti r > — oo. 
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Let (pj)clZM be a minimizing sequence. By (|5.1|l . J^(pj) is bounded, and 
by (\E'2) and the fact that / Pj = M, the minimizing sequence is bounded in 
L 1+1 / n {R 3 ). Hence we can — after extracting a subsequence — assume that it 
converges weakly to some function po £ L 1+1 /™(R 3 ), i.e., for any test function 
aeL 1+n (R 3 ), J Pj a^ J p a, cf. J, Sec. 2.18]. 

As pointed out above, the main difficulty is to prove that the induced fields 
converge strongly in L 2 ; such a result is referred to as a compactness result. It 
is true if the sequence (pj) remains concentrated: 
Step 2: Concentration implies compactness. Assume that 

lim / pj = 0, (5.2) 



J^OO 



\x\>R 



for some R>0, i.e., asymptotically the mass remains within the ball Br. We 
claim that under this assumption VU Pj — > VU Po strongly in L? . 

By weak convergence, po > almost everywhere — if po were strictly negative 
on some set S of positive, finite measure the test function a = ls would yield a 
contradiction. Moreover, Eqn. H5.2JI shows that po vanishes outside the ball Br. 
Again by weak convergence, J po = M. The sequence <jj := pj — po converges 
weakly to in L 1+1 ' n , J \oj\ < M, and O holds for \aj\ as well. We need to 
show that VU„ j — > strongly in L 2 which is equivalent to 

Ij:= ff^M dydx ^ (5 . 3) 

J J \x-y\ 

We choose R> such that Eqn. I|5.2(l applies. For e > we split the domain of 
integration into the three sets defined by 

\x-y\<e, 

\ x ~ y\ > e and (|a;| > R or \y\ > R), 
\ x ~~y\ > e an d \ x \ < R an< i \y\ < R, 

and denote the corresponding contributions to Ij by Ij t i,Ij£,Ij,3. Since 2n/(n + 
1) + 2/(n + 1) = 2, Young's inequality 7J Thm. 4.2] implies that 

l4l|<qk i ||? +1/ J|l B jr 1 |l(n + l)/2<C'n r( 3 -)/ 2 drj -0, 
for e— >0, uniformly in j. Clearly, 

\Ij.2\<-M I \<Tj(x)\dx-+Q 

6 ./|a:|>.R 

as j — > oo, for any fixed e>0. Finally by Holder's inequality, 



o~j{x)hj{x)dx 



< 



\ a j\\l+l/n\\hj\\l+n <C\\hj\\ 1+n , 
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where in a pointwise sense, 



hj(x):=l BR (x) [ l BR (y) , 1 , (Tj(y)dy->-0 

J\x~y\>e \ X ~V\ 

due to the weak convergence of Oj and the fact that the test function against 
which cTj is integrated here is in L 1+n . But since \hj\<M/e uniformly in j 
Lebesgue's dominated convergence theorem implies that hj—>0 in L 1+n , and 
the proof of l|5.3|) is complete. 

I wish to thank my student M. Hadzic for the above rather direct compact- 
ness argument |15| . 

The next two steps aim to show that minimizing sequences remain concen- 
trated and do not split into far apart pieces or spread out uniformly in space: 
Step 3: Behavior under rescaling. For p£TZm and a,b>0 we define p{x):= 
ap(bx). Then 



J pdx = ab~ 3 J pdx, E pot {p) = a 2 b- 5 E pot (p), J ^>{p) = b- 3 J tf(ap) 



dx. 



First we fix a bounded and compactly supported function p € TZm and choose 
a = b 3 so that p^TZm as well. By ($3) and since 3/n'>l, 

H r (p) = b- 3 J q{b 3 p)dx + bE pot (p)<Cb 3/n ' +bE pot (p) < 0, 

for b sufficiently small, and hence /im<0 for any M>0, i.e., a minimizcr — if 
one exists — is going to be a bound state. 

Next we fix two masses < M < M. If we take a = 1 and b = (M j~M) 1 / 3 > 1 
then for p£TZm and p£lZjj rescaled with these parameters, 

H r (p) = b- 3 J y(p)dx + b- 5 E pot (p) 

> b- 5 (J y(p)dx + E pot (p^j = (M/Mf /3 Hr(p). 

Since for the present choice of a and b the map p t— > p is one-to-one and onto 
between IZ^i and IZjj this estimate gives the following relation between the 
infima of our functional for different mass constraints: 

h w > (M/M) 5/3 h M , 0<M<M. (5.4) 

Step 4 : Spherically symmetric minimizing sequences remain concentrated. In 
this step we prove Eqn. I|5.2|l . but to make matters easier we consider for a 
moment spherically symmetric functions pGTZm, i-e., p(x) = p(\x\). For any 
radius R > we split p into the piece supported in the ball Br and the rest, i.e., 

p = Pi+ P2, Pi{x) = for |x| > R, P2{x) = for \x\ < R. 
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Clearly, 



n r (p)=n r ( Pl )+n r ( P2 )- f p -M£^L dxdy 

J x-y 



Due to spherical symmetry the potential energy of the interaction between the 
two pieces can be estimated as 



' Pi(x)p 2 (y) 
\x-y\ 



dxdy 



U Pl p 2 dx 



Air I — / pi(s)s 2 dsp2(r)r 2 dr< 
Jr r J 



(M — m) m 
R : 



where m = J p 2 is the mass outside the radius R which we want to make small 
along the minimizing sequence. We define 



Rn 



3M 2 



>0 



M 



and use the scaling estimate (|5.4H together with the fact that Km < and £ 5 / 3 - 
(1 - C) 5/3 < 1 - f £ (1 - for < £ < 1 to conclude that 



Hr(p) > 



> 



(M — m)m 
R 

( M — m\ 5 ^ 3 / m \ 5 / 3 



V M 

> h M - 



Km 



[M — m) m 
R 



1 1 

Rq R 



(M — m)m. 



(5.5) 



We claim that, if R> Rq, then for any spherically symmetric minimizing se- 
quence (pj) CUm of TC r , Eqn. (|5.2(l holds. Assume this assertion were false so 
that up to a subsequence, 



lim 

i^°°J\*\ 



J 

J\x\>R 



Pj 



--m>0. 



Choose Rj > R such that 



J\x\>R? z J\i 



\x\>R 



By |53|, 

H r (pj) > h M 



1 1 

Ro Rj 



(M — rrij ) rrij > hu + 



1 1 

Rq R 



and letting j — > oo leads to a contradiction. 

Step 5: Removing the symmetry assumption. The restriction to spherical sym- 
metry would mean that stability would only hold against spherically symmetric 
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perturbations. Fortunately, the restriction can be removed using a general re- 
sult due to Burchard and Guo. To explain it we define for a given function 
p£L+(M. 3 ) its spherically symmetric decreasing rearrangement p* as the unique 
spherically symmetric, decreasing function with the property that for every 
r>0 the sup-level-sets {xeR 3 \p(x)>T} and {ieR 3 |p*(i) >t} have the same 
volume; the latter set is of course then a ball about the origin whose radius 
is determined by the volume of the former. The important point is that for 
any monotone function "J" the integral J ^(p)dx does not change under such a 
rearrangement, while the potential energy can only decrease, and it does not 
decrease if and only if p is already spherically symmetric (with respect to some 
center of symmetry) and decreasing. These facts can be found in (7J Ch. 3]. In 
particular, a minimizcr must a posteriori be spherically symmetric. 

Now let (pj) CIZm be a not necessarily spherically symmetric minimizing 
sequence. Obviously, the sequence of spherically symmetric decreasing rear- 
rangements (pj) is again minimizing. Hence by the previous steps, up to a 
subsequence (p*) converges weakly to a minimizer po = /°o an d 

VI7pj -> VC/ in L 2 , hence / * (p* ) -> / * (p ) . 



Moreover, 



-> Hr(po)- I y{po)=E pot ( Pa ). 



In this situation the result of Burchard and Guo ^3 Thm. 1] says that up to 
translations in space 

VU Pj — > VU in L 2 (5.6) 

The proof of this general result is by no means easy, and it is possible to 
obtain stability against general perturbations without resorting to it, cf. [B1IH1 
I14j . However, since this general result may be useful for other problems of this 
nature we wanted to mention and exploit it here. 

Step 6: Proof of Theorem^ Given a minimizing sequence (pj) we already know 
that up to a subsequence it converges weakly in L 1+1 / n to a non-negative limit 
po of mass M. The functional p\— > j^(p)dx is convex by Assumption (^1), so 
by Mazur's Lemma Thm. 2.13] and Fatou's Lemma Lemma 1.7] 

/ *(/9 )(ia;<limsup / ^(pj)dx } (5.7) 

J j^oo J 

in particular, po £1Zm- Together with Eqn. H5.6(l this implies that 

Ti-r (po ) < limsup7i r (pi) = h M 

i — >oo 

so that po is a minimizer of Ti. r , and the proof of Theorem is complete. 
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6 Related results, concluding remarks, and 
open problems 

In this last section we want to touch upon a number of questions which may 
come to mind. 

The threshold k — 3/2. First we want to ask what happens if we choose the 
exponent k in the assumptions ($2), ($3) larger than 3/2. This question is 
answered by the following observations: For the Vlasov-Poisson system the 
Casimir functional C is preserved. Hence we can pursue an alternative approach, 
namely to minimize the energy 

"H = Ekin + Epat 

under the mass-Casimir constraint 

J I fdvdx + C(f) = M. 

As was shown in jHUHj this works, provided $ € C^QOjOoQ with $(0) =0 = $'(0) 
satisfies the assumptions ($1) (strict convexity) and 

($2') $(/) > C/ 1+1 / fe for / > large, with < k < 7/2; 

($3) is not needed any longer. However, reduction, which combined the kinetic 
energy and the Casimir functional into the new functional J ^S(p), can no longer 
work, because the former two functionals now appear in different places in the 
variational problem. Moreover for polytropes one can show that the energy- 
Casimir functional changes its sign from negative to positive at n = 3 i.e. k — 3/2. 
If a perturbation leads to an initial datum with positive energy 7i(/(0)) > then 

sup{|x| | (x,v) E supp/(t)} > Ct, i— >-oo. 

The analogous result holds for the Euler-Poisson system, except that for a min- 
imizer the energy in the Euler-Poisson picture equals the energy-Casimir func- 
tional in the kinetic Vlasov-Poisson picture. Hence the fact that He (for the 
polytropes) changes sign at fc = 3/2 does not signify any stability change on the 
kinetic level, but it does destroy stability on the fluid level. 
Are all relevant isotropic models covered? To answer this question it is useful to 
translate the conditions on the function $ determining the Casimir functional 
into conditions on the dependence fo(x,v) = (j)(E — E) of the resulting steady 
state on the particle energy. Since ($') _1 = 0, $ is strictly convex as required 
iff strictly increases; if <p(n) < Cn k for rj large and <f>(rj) > Cr/ k for i] > small 
then $ satisfies the growth conditions ($2) and ($3), provided 0<fc,fc'<3/2. 
In order to satisfy the alternative growth condition ($2') we only need to require 
0(^)<C?? fc for rj large with 0</c<7/2. This means that all those polytropic 
steady states which are compactly supported and decrease as a function of the 
local energy are covered by our results — either by the approach via reduction 
with the bonus of the stability of the corresponding fluid model or, when this is 
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no longer possible, by the alternative approach mentioned above. The alterna- 
tive approach also covers the limiting case f (x,v) = (E -E) 7 / 2 , the so-called 
Plummer sphere, cf. |SJ. The minimization property of the latter is also inves- 
tigated in ^2] i without deducing stability. 

An important example from astrophysics which is not yet covered is King's 
model where fo(x,v) = exp(E — E) — 1 on its support. It leads to a Casimir 
function $ growing like / In/ which is to slow. Possible extensions of the method 
to include King's model are currently under investigation. 
Does the method apply to non-isotropic steady states? In the presence of spher- 
ical or axial symmetry (with respect to the X3-axis) angular momentum quan- 
tities like 

L := \x x v\ 2 or L3 := X1V2 — X2V1 = (xx 77)3 

are conserved along particle trajectories. If we let the function $ depend also 
on L or L3 the corresponding variational problem has a solution which is a 
steady state depending on the additional invariant. The dependence on L was 
treated in ^3 ^] [201 while axially symmetric steady states depending on L3 
were treated in |21| . The main assumption again is the strict monotonicity of the 
dependence on the local energy. There is however one important difference to 
the isotropic case: Since a Casimir where $ also depends on angular momentum 
is preserved by time dependent solutions only if they have the corresponding 
symmetry we obtain stability only with respect to either spherically symmetric 
or axially symmetric perturbations. The question of stability of non-isotropic 
steady states against non-symmetric perturbations is under investigation. 

The method has also been applied to disk-like steady states, cf. [221 ■ Here 
we had to assume that the perturbations live only in the plane of the disk. An 
extension of these results in the spirit of the later developments reported in the 
present notes or in [S] is under investigation. Stability against perturbations 
perpendicular to the disk is another challenging open problem in this area. 
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